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The tetragonal states produced by isotropic pseudomorphic epitaxial strain in the (001) plane on 
a tetragonal phase of a crystal are calculated for V, Ti, Rb, Li, K, Sr from first-principles electronic 
theory. It is shown that each metal has two tetragonal phases corresponding to minima of the 
total energy with respect to tetragonal deformations, hence are equilibrium phases, and that the 
equilibrium phases are separated by a region of inherent instability. The equilibrium phase for any 
strained tetragonal state can thus be uniquely identified. Lattice constants and relative energies of 
the two phases and the saddle point between them are tabulated, as well as the tetragonal elastic 
constants of each phase. 



I. INTRODUCTION 

Metal crystals in the body-centered tetragonal struc- 
ture are known to have two total-energy minima as 
functions of tetragonal lattice constants a and c 
The minima are appropriately called equilibrium metal- 
lic phases, since they persist without applied stress and 
are stable under small tetragonal deformations. Previ- 
ous work by the authors Q defined and discussed epi- 
taxial Bain paths (EBP), which are sequences of tetrago- 
nal states that include the equilibrium tetragonal phases. 
These paths give the strained tetragonal states pro- 
duced by isotropic epitaxial strain, i.e., isotropic two- 
dimensional biaxial or in-plane strain, imposed on the 
(001) planes of the equilibrium phases. In the previous 
work the EBP for V, Co, and Cu were found using first- 
principles total-energy calculations and the EBP were 
compared with the paths produced by uniaxial stress on 
the phases, which were called in that work uniaxial Bain 
paths. 

The EBP were shown to be directly useful in interpret- 
ing the bulk structure of epitaxial films determined, for 
example, by quantitative low-energy electron diffraction 
(QLEED). Comparison of the measured structure of a 
film with the states on the EBP identifies the phase from 
which the film is produced by the epitaxial strain. Thus 
a film of Co on Cu(001) was shown to be strained fee Co, 
but a film of Co on Fe(001) was shown to be strained 
body-centered tetragonal (bet) Co, a metastable phase 
of Co, whereas bec Co was shown to be unstable. 

In the tetragonal plane, whose coordinates are tetrag- 
onal lattice constants, the EBP is a continuous path that 
passes through the two phase points at the energy min- 
ima and through the saddle point of energy between the 
two minima. It was shown that between the minima 
a segment of the EBP exists which includes the saddle 
point, but not the minima, that consists of inherently 
unstable states. More generally, strained states of each 
phase, not just those on the EBP, were shown to be sepa- 
rated by a region of inherently unstable tetragonal states, 
so that an observed strained tetragonal state has a con- 



nection to just one equilibrium phase through stable but 
constrained states. 

The present work gives the EBP for six metals based 
on the published first-principles calculations for tetrago- 
nal structures by Sliwko, Mohn, Schwarz and Blaha [||. 
The calculation procedures for finding the EBP, the con- 
tours of constant energy and the unstable region are de- 
scribed in Sec. II. The results are described in Sec. Ill with 
two tables and four figures. Section IV discusses why the 
EBP is useful, the significance of the unstable region, and 
notes defects and generalizations of these calculations. 

II. CALCULATION PROCEDURES 

The calculations of total energy E. as a function of 
the tetragonal lattice parameters a, the side of the square 
cross section, and c, the height of the unit cell, used the 
power-series expansions given in Ref . , whereas the cal- 
culations in Ref. || used the full-potential APW program 
Wien95 directly. The power-series expansions were fitted 
in Ref. Q to extensive first-principles calculations in the 
local-density approximation (LDA) with Wien95. The 
expansions give E within specified ranges of c/a and vol- 
ume per atom V ~ ca 2 /2 that include the minima and 
saddle point and have the form 

n m 

J? = ^^i y ( C / fl )'P. (1) 

i=0 j=0 

The coefficients Aij are tabulated in Ref. || to eight sig- 
nificant figures and are available in electronic form from 
the authors. In Eq.(|l|) n is 5 or 6, m is 3 and E is ob- 
tained in the specified ranges to an accuracy stated to 
be better than 0.01 mRy ||. Some comments on this 
stated accuracy are made in Sec. IV. A useful feature of 
the formula (1) is that analytical formulas for the first 
and second derivatives of E may be readily derived. 

The EBP for each metal is found by calculating E, 
(dE/dc) a , and (dE/da) c as a function of c at constant 
a, and locating the c value for which E has a minimum 
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or, equivalently, locating the zero of (dE/dc) . The min- 
imum corresponds to the epitaxial film condition of zero 
normal stress on the (001) surface. As a ranges over the 
structures between the phase points, the EBP is traced 
by these minima of E at each a. Since E and (dE/dc) a 
are evaluated easily, a dense grid of c values permits sim- 
ple interpolation to four significant figures for c,E, and 
(dE/da) c at the minimum of E for any a. 

The two phase points and the saddle point correspond 
to stationary points for E, hence are located by interpo- 
lating on the EBP itself to find points where (dE/da) c 
vanishes along with (dE/dc) a . Two of the three station- 
ary points are always cubic points, since at a cubic point 
if (dE/dc) a = then also {dE/da) c = 0. The contours 
of constant E are similarly found by interpolating the 
desired E in a tabulation of E(c) at values of a over a 
range of a that covers the desired contour. 

Tetragonal elastic constants at any a and c may be 
defined by 
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The Cij differ from the usual elastic stiffness coefficients 
Cij because the cy correspond to tetragonal deformations 
which maintain the square symmetry a\ = a?, || . To sep- 
arate en from C12 requires breaking tetragonal symmetry, 
but is not possible if E is known only from the power se- 
ries (1). However when the phase has cubic symmetry 
en = C33 and C12 can be evaluated. In fact :ci2 can be 
evaluated in two ways, i.e., from en and C33 on the one 
hand, and from 613 on the other hand. The correspon- 
dence of the two values is then a test of the accuracy of 
the power series representation of E as will be noted in 
Sec.IV. 

A strained tetragonal state will in general be main- 
tained by applied in-plane and out-of-plane stresses de- 
termined by the derivatives (dE/da) c and (dE/dc) a . 
However stability depends also on a condition on the sec- 
ond derivatives of E which states that the second-order 
differential of E is always positive, i.e., that 
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is greater than zero for all deformations 5a and Sc. Oth- 
erwise the structure would have a tetragonal deformation 
that lowers the energy, so that the structure cannot be 
maintained by applied stresses. The condition 5 2 E > 
is then a condition on the cy , namely, 



The lines along which D = can be calculated readily by 
finding the c/a at which D = for the function E(c/a) 
at constant V and using a range of V to follow the line; 
the analytical power series for the second derivatives of 
E. obtained from (1) are convenient for the calculation. 

The slope of the EBP at the phase points can be ex- 
pressed directly in terms of the elastic constants of each 
tetragonal or cubic phase, volume Vq = c^a^/2, 



d(V/V ) 
d(c/a) 



EBP 



Qq (2c 33 - ci 3 ) _ 2a (c 33 - ci 3 ) 
Co (C33 + C13) c (c 33 + 2ci 3 )' 

(5) 



If the phase point has cubic symmetry, Eq.(gj) simplifies, 
since then C\\ = C33, C12 = C13 and 
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where v is the Poisson ratio of the cubic phase, v = 
ci2/(cn + C12). In the figures Vq is chosen as the vol- 
ume of the equilibrium phase of lower E, hence Eq.(|^) 
or Eq.(^|) applied to the other equilibrium phase has a 
factor of the ratio of the V of the other phase to Vq on 
the right sides of Eq.(||) and Eq.(^j|) to get the slope of 
the plotted EBP. 

Equations (||) and (|^) relate the linear elastic approxi- 
mation to the EBP of a phase directly to the elastic con- 
stants for a tetragonal noncubic phase or a cubic phase 
respectively. If the elastic constants of the phase are 
known from experiment or theory, the equations give the 
linear approximation to the EBP. Then the measured 
bulk structure of a strained epitaxial film can be com- 
pared with the linear EBP to identify the equilibrium 
phase of that film. This identification is especially in- 
teresting for noncubic tetragonal phases, which are pre- 
dicted to exist for all transition elements but are 
always metastable. Hence they cannot be made macro- 
scopically, but may be stabilized by epitaxy, as was done 
i|. Comparison of Eq.(||) with 
) defines an effective Poisson ra- 



in the case of bet Co 
Eq. (^|) shows that Eq. (|_ 
tio for tetragonal phases 
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III. RESULTS 



The results of calculation with the procedures and for- 
mulas of Sec. II are given in four figures and two tables. 
Figures 1 to 3 plot the EBP of V, Ti, and Sr along with 
contours of constant E on the c/a — V/Vo tetragonal 
plane, where Vq is the volume per atom of the more sta- 
ble phase point; the positions of the two phase points and 
the saddle point are marked. The composite Fig. 4 plots 
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the EBP of Rb, Li, and K without the contour lines. The 
corresponding coordinates and the energy at each point 
referred to a zero at the more stable phase are given in 
Table I. Plotted in Figs.l, 2 and 3 are five contours of 
constant energy, i.e., two contours at SE above the two 
minima, the contours through the saddle point, and the 
contours SE above and below the saddle-point energy. 
The values chosen for SE depend on the energy scale for 
each metal. The unstable region where D < 0, which 
includes the saddle point, is the region between the two 
lines of long dashes, and is shown in all four figures. Ta- 
ble I also gives the stationary points of E for tetragonal 
Rb computed by Milstein, Marschall and Fang from 
an empirical potential fitted to experiment. 
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FIG. 1. EBP for V on the c/a - V/V plane (full line), 
phase point 1 bcc (full circle), volume Vb = 84.8bohr 3 , 
E = 0; saddle point fee (marked x) E = 20.8 mRy; phase 
point 2 bet (full circle), E = 10.0 mRy. Energy contours 
are drawn (short dashes) at SE — 1 mRy above the min- 
ima, through the saddle point and ±SE from the saddle 
point. The unstable region is between the lines of long 
dashes and includes the saddle point. The coordinates 
for the phase and saddle points are in Table I. 

The tetragonal elastic stiffness constants en , C13, C33 
are given in Table II at the phase points and also the 
elastic stiffness constants found in Ref. || for the cu- 
bic phases of each metal. When the phase is cubic, the 
usual tetragonal elastic constants en and C\i (found in 
two ways) are given; note that then en = C33. By sym- 
metry, E is stationary at the cubic points (c/a = 1, \/2) 
on the EBP, i.e., dE/da = dE/dc = at cubic points. 
However the cubic points can be saddle points as well 
as minima. The possible configurations have been clas- 
sified in Ref. Q in three cases, i.e., Case 1: minimum 
at bcc (c/a = 1), saddle point at 1 < c/a < \/2, and 
minimum at fee (c/a = \/2); Case 2: minimum at bet 
(c/a < 1), saddle point at bcc (c/a = 1), minimum at fee 
(c/a = V2); Case 3: minimum at bcc (c/a = 1), saddle 



point at fee (c/a = V2), minimum at bet (c/a > \/2). 
Then Rb, K, Li, Sr are Case 1, Ti is Case 2 and V is Case 
3. 
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FIG. 2. EBP for Ti, stationary points, contour lines 
and unstable region marked as in Fig. 1. Phase point 1 
is bet at E = 1.59 mRy; saddle point is bcc at E = 2.59 
mRy; phase point 2 is fee at E = mRy, Vo — 108.1 
bohr 3 ; SE = 0.2 mRy. The measured strained bulk struc- 
ture of epitaxial film on Al(001) is marked by the open 
square with error line 




FIG. 3. EBP for Sr, stationary points, contour lines 
and unstable region marked as in Fig, 1, Phase point 1 
is bcc at E = mRy; saddle point is bet at E = 0.383 
mRy; phase point 2 is fee at E = 0.190 mRy, V a = 319.6 
bohr 3 , SE = 0.07 mRy. 

The contour lines of constant E 7 which are vertically- 
oriented ellipses near each phase point, also appear in 
Ref. H plotted on the c/a — V plane for each metal. The 
EBP and contour lines for V are also in Ref. Q, where 
they are plotted on the a/ao — V/Vq plane. On this plane 
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the contours are tilted and the bcc and fee positions on 
the EBP are reversed. The contours on the c/a—V plane 
in terms of the deviations of c/a and V from the values 
at the phase points are given by 
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FIG. 4. EBP for Rb, Li, K, stationary points, and un- 
stable region marked as in Fig. 1. Phase point 1 for Rb is 
bcc at E = 0.135 mRy; saddle point is bet at E = 0.157 
mRy; phase point 2 is fee at E = mRy, Vb = 518.0 
bohr 3 . Phase point 1 for Li is bcc at E = 0.144 mRy; 
saddle point is bet at E = 0.161 mRy; phase point 2 is 
fee at E = mRy, V = 127.7 bohr 3 . Phase point 1 
for K is bcc at E = 0.013 mRy; saddle point is bet at 
E = 0.089 mRy, V = 430.4 bohr 3 . 

For cubic symmetry ci vanishes; even for the bet 
phases of Ti and V, ci is much smaller than c\ and c 3 , 
so the contours are still nearly vertical. Formula (7) for 
the case of cubic symmetry with coefficients in terms of 
the Cjj is given in Ref. . 



IV. DISCUSSION 

The principal result found here is the EBP between the 
phase points, which shows the effects of isotropic epitax- 
ial strain on equilibrium phases. The EBP provides a 



basic elastic response of a material in tetragonal struc- 
ture to a particular strain, one which is available experi- 
mentally, including the interesting case of negative strain 
in the plane. The EBP are found here from first princi- 
ples in a well-defined approximation with errors of known 
magnitude, and include nonlinear effects. These EBP can 
then be compared directly with the strained bulk struc- 
tures determined by QLEED for epitaxial films. This 
comparison is illustrated for strained epitaxial Ti films 
on Al(001) Q, where the QLEED point and its error 
limits (from the uncertainty in the bulk value of c) are 
shown in Fig. 2 to agree well with the calculated EBP. 
The Vb used in evaluating V/Vq for the QLEED point is 
the measured hep Ti value, V hcp — 119.2bohr 3 which is 
close to the fee value. 

The presence on the EBP of an inherently unstable sec- 
tion separating strained fee Ti from strained bet Ti is an 
important result of the theory. Since the QLEED point 
within its error limit lies on the fee part of the EBP, the 
epitaxial film must be strained fee Ti. This result is par- 
ticularly interesting because fee Ti does not appear on the 
usual pressure - temperature phase diagram. Note that 
the theory gives directly the possible states of epitaxially- 
strained fee Ti, including any nonlinear elastic behavior 
of the crystalline phases. This comparison of measured 
structure with the theoretical EBP replaces the previous 
analysis, which assumed constant elastic stiffness coeffi- 
cients and attempted to estimate the elastic coefficients 
of the cubic phases |Q . The use of linear elastic relations 
for identification of the equilibrium phase is compared to 
the use of the EBP for that identification in the case of 
Co in Ref. 0], where tetragonal states and the EBP are 
plotted on the a — d plane, where d is the layer spacing. 

These first-principles calculations, which do not use 
empirical information, have an advantage over empirical 
potentials fitted to measurements, since these calcula- 
tions are as good for highly strained or even unstable 
states, which are not accessible to measurement, as they 
are for slightly strained equilibrium states. Thus they 
provide a test of calculations based on empirical poten- 
tials, such as the potential used for Rb in Ref. (t]]. Com- 
parison in Table I shows that the c/a values of the saddle 
point differ by 10%, that the minimum energy in Ref. 
is at the bcc structure, rather than the fee structure 
found here, and that the energy separation of the equilib- 
rium phases in Ref. has both a sign and a magnitude 
different (smaller) from what is calculated here. 

The power series representation of the tetragonal en- 
ergies shows some defects, e.g., the cubic points deviate 
from c/a = 1 or \/2. in Table I, the two values of C12 
do not agree well in some cases. These defects are in the 
representation of the results of LDA calculations. The 
most serious defect is in the saddle point of Ti, which 
is bcc, but the power series finds the saddle point at 
c/a = 1.05 rather than 1. A recalculation of the energy 
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E of Ti directly with WIEN95 finds that the power series 
has missed an asymmetry in E around the saddle point 
which shifts the position of the maximum. The direct 
calculation finds the maximum at c/a = 1 as it should 
be; it also verifies the minimum of fee Ti at c/a = 1.40. 
The deviation from symmetry at fee Ti and all other cu- 
bic phases is thus no more than 1%, except for fee Rb, 



where the deviation is 1.7%. In comparison with experi- 
ment all the calculated elastic stiffness coefficients are too 
large by at least 10% and the volumes per atom are too 
small by 5 to 10%. These discrepancies from experiment 
are defects of the assumptions of the band calculations, 
i.e., of the LDA with semi-relativistic corrections. 



TABLE I. Tetragonal states stationary in energy. The parameter are (c/a)i,(V/Vo)i, Ei, i — 1,2,3, in the tetragonal plane 
for states of stationary E; the stable phase point is i — 1 or 3 and the metastable phase point is then i = 3 or 1, respectively; the 
saddle point is i = 2. Note that the tetragonal lattice parameters can be found from a = [2V (V/V )/(c/a)]^ 3 and c = a(c/a). 
Vo is the volume per atom of the more stable phase point in bohr 3 and Ei is the energy in mRy with respect to the energy of 
the more stable phase point. 



Metal 


Vo 


(c/a)! 


[V/Vo)i 


E x 


(c/a) 2 


{V/Vo) 2 


E 2 


(c/a) 3 


(V/V ) 3 


E 3 


V 


84.9 


1.00 


1.00 


0.00 


1.43 


1.04 


20.8 


1.78 


1.02 


10.0 


Ti 


107.9 


0.85 


0.99 


1.59 


1.05 


0.99 


2.59 


1.40 


1.00 


0.0 


Rb a 


518.0 


0.99 


0.999 


0.135 


1.10 


1.001 


0.157 


1.39 


1.00 


0.0 


Rb 6 




1.00 


0.997 


0.042 


1.22 


0.999 


0.041 


1.41 


1.00 


0.0 


K 


430.4 


1.00 


0.996 


0.013 


1.19 


1.000 


0.089 


1.42 


1.00 


0.0 


Li 


127.7 


0.99 


1.003 


0.144 


1.10 


1.003 


0.161 


1.41 


1.00 


0.0 


Sr 


319.6 


1.00 


1.000 


0.000 


1.21 


1.009 


0.383 


1.42 


1.014 


0.109 



a This work. 
b From Ref. §i 



TABLE II. Elastic constants of stable and metastable phases: The tetragonal elastic stiffness constants Cn, C13, C33, in 
Mbar for phase 1 at stationary point 1 of Table I and for phase 2 at stationary point 3 of Table I, Cn, C13, C33 are the usual 
elastic constants for tetragonal structures. For cubic phases en = C33 = C33, hence from Eq.(2), C12 is found in two ways: 
c i2 = (cn/2) - C33 and c$ = cis/2. Note that 1 mRy/bohr 3 = 0.14711 Mbar. 



T7T1 / TO TO / TO TO~~ 

Metal c/a en C13 C33 c\ 2 c\ 2 c/a C11 C13 C33 c\ 2 c\ 2 

1.00 10.26 3.27 3.60 1.53 1.63 1.78 11.47 2.08 4.57 

V 6 2.88 1.36 

Ti a 0.85 5.05 2.18 1.45 1.40 5.52 2.38 1.59 1.17 1.19 

Ti 6 1.48 1.21 

Rb a 0.99 0.160 0.074 0.043 0.037 0.037 1.39 0.163 0.075 0.048 0.034 0.038 

Rb 6 0.045 0.038 0.046 0.040 

K a 1.00 0.220 0.100 0.059 0.051 0.048 1.42 0.214 0.097 0.062 0.046 0.049 

K b 0.061 0.049 0.057 0.050 

Li a 0.99 0.633 0.296 0.165 0.151 0.148 1.41 0.636 0.280 0.178 0.140 0.140 

Li 6 0.170 0.147 0.169 0.143 

Sr a 1.00 0.678 0.283 0.194 0.145 0.141 1.42 0.698 0.295 0.202 0.147 0.147 

Sr 6 0.206 0.136 0.189 0.153 



a This work. 
6 From Ref. §. 



Despite the deviations from experiment, which may 
be reduced in subsequent calculations by more accurate 
formulation of the electronic structure equations, these 
results are of immediate practical value in interpreting 
measured film structures, and of conceptual value in pro- 
viding a sharp distinction between, for example, a tetrag- 
onally strained bee phase and a tetragonally strained fee 
phase. A generalization of the tetragonal results to other 



structures suggests a new type ofphase diagram in which 
the various equilibrium phases, stable and metastable, 
are points in a parameter space which has structural pa- 
rameters as coordinates. The present results suggest that 
each phase point is surrounded by a region of strained 
states and the regions are embedded in and separated 
by a continuous matrix of inherently unstable states that 
cannot be stabilized by applied stresses. Such a general- 
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ization for structures with considerable symmetry, such 
as the tetragonal structure, seems calculable by present 
codes. 
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